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We establish a set of nonequilibrium quantum phase transitions in the Lipkin-Meshkov-Glick 
model driven under monochromatic nonadiabatic modulation of the transverse field. We show that 
the external driving induces a rich phase diagram that characterizes the multistability in the system. 
Interestingly, the number of stable configurations can be tuned by increasing the amplitude of the 
driving field. Furthermore, by studying the quantum evolution , we demonstrate that the system 
exhibits a set of quantum phases that correspond to dynamically stabilized states. 
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A quantum phase transition (QPT) is a drastic change 
of state of a many-body system which occurs at zero 
temperature In contrast to thermodynamic phase 
transitions — which are driven by thermal fluctuations, 
QPTs are entirely driven by quantum fluctuations. This 
critical behavior gives rise to states of matter, which — 
in some cases, can be described in terms of symmetry- 
broken phases. Most of the research in the past has 
been focused on QPTs on equilibrium systems, whereas 
relatively little is known about the critical behavior 
in nonequilibrium conditions. The quantum control of 
many-body systems by a driving field has attracted con- 
siderable interest both theoretical and experimental with 
workers from very different communities beginning to 
look at driven models @-GJ]. The possibility to manip- 
ulate the quantum state of a system by means of a clas- 
sical external control allows one to explore novel states 
of matter and effective interactions which are absent in 
equilibrium [Tl| - [T3 |. In this context, a recent line of re- 
search are excited-state phase transitions, e.g. in mean- 
field type QPTs such as the Lipkin-Meshkov-Glick [l5[ 
or the Dicke superradiance [l|| model. Furthermore, the 
connection between chaos and criticality fl7| has been 
studied in detail for the interacting boson model (l8j . 

Despite of the original motivation of the Lipkin- 
Meshkov-Glick (LMG) model as a toy model to test 
approximation methods in manybody physics [l9| . cur- 
rently it constitutes an active field of research and a nat- 
ural scenario to study the relation between QPTs and 
spin squeezing [2(| and quantum Fisher information as 
a resource for high-precision quantum estimation 21 1. 
Furthermore, rather recently, experimental realizations 
of the LMG model in optical cavity QED [H, [H] and 
circuit QED [24j have been suggested. 



In this paper we study a driven version of the LMG 
model in which we assume a time-dependent inter- 
particle interaction. In comparison with previous works 
which considered dynamical properties of the LMG 
model under an adiabatic change of parameters of the 
system across the quantum critical point (2||, the ef- 



fect of fast and slow quenches of the transverse field [26| , 
and a periodically-driven uniaxial LMG model [q. |27I. |28| . 
we address here the fundamental issue of driving-induced 
QPTs in the LMG model. 

The development of novel experimental techniques to 
externally control many-body systems has been exten- 
sively increased in the last years. For example, the kicked 
top, has been realised experimentally in an ensemble of 
laser-cooled Cs atoms [29(. Furthermore, an optical re- 
alization of the uniaxial driven LMG in photonic lattices 
[3Cl | and by using superconducting charge qubits con- 
nected in parallel to a common superconductor induc- 
tance [3l| have been suggested. 

We show that when the driving is near resonance with 
the excitation energy of the undriven system the nature 
of criticality changes dramatically. Additionally, we find 
a phase diagram replete with a host of macroscopically 
distinct metastable phases with no analogues in equilib- 
rium. This gives rise to a novel route of experimental 
studies exploring the characteristics of criticality under 
nonequilibrium conditions. 

A more specific outline of our paper is as follows. In 
section fl] we describe the fundamentals on the QPT in 
the undriven LMG model and construct a bosonized ef- 
fective Hamiltonian which allows us to describe the sta- 
bility properties of the symmetric phase. In section ILT1 we 
find an effective Hamiltonian by using the rotating wave 
approximation (RWA) and describe the quantum evolu- 
tion of the observables to understand the characteristics 
of the novel nonequilibrium metastable phases. Finally, 
conclusive remarks arc given in section IPLT1 



I. QUANTUM RESONANCES IN THE LMG 
MODEL. 



The periodically-driven Lipkin-Meshkov-Glick model 
describes the dynamics of N interacting two-level sys- 
tems in a transverse local field 



H(t) = -hJ z - 1 (-f{t)J a 
,(*) 
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where J a — h J2i=i °a denote collective angular mo- 

(i) 

mentum operators and <r a are the Pauli matrices with 
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a G {x,y, z}. These operators satisfy the SU (2) al- 
gebra [J a ,Jp] = itafi-yJ-y- In the following we shall 
consider ft < 0, and a monochromatic modulation of 
the interparticle interaction with a static contribution: 
7 x (t) = 7q + 7f cos fit. 

In this paper we consider Dicke states i.e., the whole anal- 
ysis is reduced to the Hilbcrt subspacc characterized by a 
maximal total angular momentum j = N/2. Associated 
with the Hamiltonian Eq. (|T|) is a conserved parity 

ft = exp(i7r(J,+i)), (2) 

such that [H(t), II] = 0. Our aim in this paper is to study 
the new aspects of criticality under the effect of driving. 
In this section we provide the basics on the formalism 
used to describe the effective Hamiltonian for the sym- 
metric phase of the LMG model. In particular, we find a 
resonance condition related to m-photon processes under 
the effect of driving. 

A. The QPT in the undriven LMG model 

In the case of the undriven LMG model (jf = 0), 
an analytic study of the ground-state energy surface in 
the thermodynamic limit TV — > oo leads to a phase di- 
agram in the (tq , 7 y )-plane, which is divided into four 
regions depending on the geometry of the surface. These 
sections are distinguished from each other by the num- 
ber of minima, maxima and saddle points, which are re- 
lated to non-analyticities of the density of states [l5[. 
As a consequence of the symmetry of the LMG Hamilto- 
nian, it is sufficient to consider the region with |tq | < 7 y . 
This analysis shows that in region |tq| < j y < —ft the 
ground-state energy landscape has a single global mini- 
mum, whereas in the regions with \ < —ft < j y the 
surface has two global minima. By crossing the criti- 
cal line 7 y = —ft, the single global minimum splits in 
two global minima and the system exhibits a continous 
transition from a symmetric state to a symmetry-broken 
state, i.e., a second-order QPT. 



and take the thermodynamic limit N — > oo, assuming 
b/N 0. The result is a bosonized Hamiltonian for the 
symmetric phase 

H s (t) = -htfb-- A [ 1 %t){tf+b) 2 - 1 y{b^b) 2 \+^. (4) 

By introducing the coordinate operators 




we obtain the Hamiltonian of a parametrically-drivcn 
harmonic oscillator [32j 

Hs(t) = £ + \ (V + ft 7l * (l + cosfit) f + ^, 

(7) 

where 

— VK)K) < 8 > 

is the characteristic energy scale of the system in the 
absence of driving (7^ = 0). The undriven system ex- 
hibits a critical behavior which is related to softening 
of the collective excitation spectrum, i.e., when the sys- 
tem is close to the critical point (7 s — » —ft), the sys- 
tem exhibits a gapless excitation above the ground state 
(e —¥ 0). Therefore, in the region |7q| < —ft < j y the 
symmetric phase became unstable. Interestingly, in the 
driven system the situation can be highly nontrivial as a 
consequence of mechanisms such as parametric resonance 
and parametric stabilization, which are characteristic of 
the parametrically-driven harmonic oscillator [UHII • As 
a consequence, one can tune up conveniently the param- 
eters close to resonance in order to manipulate the sta- 
bility of the system, i.e., to produce a change of phase 
for parameters which are far away from the critical point 
of the undriven system. 



B. Effective bosonized Hamiltonian for the 
symmetric phase 

We begin our analysis by investigating the stability 
of the symmetric phase under the effect of an external 
driving. To this end, we construct a symmetric phase 
effective Hamiltonian in the same way as in Ref. hj : we 
make a Holstein-Primakoff representation of the angular 
momentum algebra in terms of bosonic operators b, b' 

N 

J* = tfb- Y , (3) 

J+ = - fctfc, 

J- = - tfb b, 



C. Resonance conditions 

In the thermodynamic limit, the undriven LMG model 
is characterized by a single collective excitation with en- 
ergy e. Under the effect of an external drivin g, th e possi- 
bility of multiphotonic resonances arises [33j-|35[. In the 
scmiclassical theory of light-matter interaction, we can 
interpret a Floquet state as a light-matter quantum state 
containing a definite, though very large, number of pho- 
tons [33|. Multiple transitions between quantum states 
of system that are not directly coupled by the interaction 
can occur by means of intermediate states with a differ- 
ent number of photons present [331 134. l36| ■ In particular, 
?7i-photon transitions occur when the condition 

2e = mO, (9) 
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FIG. 1: (Color online) Phase diagram of the nonequilibrium 
QPT in the Lipkin-Meshkov-Glick model. The number of 
local minima of the quasienergy surface Eq (Q, P) in the col- 
ored zones are indicated by the labels, (a) Depicts the phase 
diagram as a function of 7 y and the driving amplitude 
for 7o/|/i| = 0.5, and (b) as a function of 70 and 71 for 
7 ! '/|/i| = 2. The dashed red (light gray) line in (a) resem- 
bles the second-order QPT in the undriven LMG model that 
occurs at 7 a /|/i| = 1, in the driven case, however, this line is 
the separatrix between the regions with even and odd num- 
ber of minima. The dashed black lines in (a) and (b) depict 
the transition of the symmetric phase (Q, P) — (0, 0) from 
a saddle point to a local maximum, and correspond to the 
contour A2 = 0. We consider the parameters Q/\h\ = 40 and 
h/\h\ = -1. 



with integer m is satisfied. For a parametric oscillator 
with fundamental frequency e, Eq. (|9]) is the usual reso- 
nance condition [32} . In Floquet theory, Eq. © implies 
the existence of a crossing between the energy levels when 
the energy spectrum of the undriven system is folded into 
the Brillouin zone [33]. In this paper we focus on the 
parameter regime 5^ m \jQ,j y <C fi, where the detuning 



5<- r 



describes how far the system is from 



resonance, and the resonance condition Eq. (jH]) reads 



mf2 



(10) 



Such resonance condition will be used in the next section 
to perform a description of the system based on an ef- 
fective time-independent Hamiltonian which is valid for 
parameters close to a multiphotonic resonance. 



II. THE ROTATING WAVE APPROXIMATION 
AND THE EFFECTIVE HAMILTONIAN 
APPROACH 

As we are interested in the asymptotic quantum dy- 
namics and the description of critical signatures, it is 
convenient to describe the dynamics of the system in a 
rotating frame. In the limit S^ m \ Jq , j y <C 0, it is possi- 
ble to neglect the fast oscillations in the rotating frame, 
and a treatment of the system based on the description 
for time-independent systems is possible via an approxi- 
mate effective Hamiltonian. Motivated by the m-photon 



resonance condition Eq. (|10p . we perform a study of 
the system based on the rotating wave approximation 
(RWA) [TH, H(| . Let us perform an unitary transforma- 
tion of Hamiltonian Eq. (JTJ) into a convenient rotating 
frame via the unitary operator 

U m {t) = exp (-i0(t) Jl) exp (—iO m (t)J z ) , (11) 

where 9(t) = ^ and 9 m (t) = sfL. In the rotat- 
ing frame the dynamics is governed by the Hamiltonian 
H m {t) = Ul(t)HUrJt), where H = H{t)-i§- t is the Flo- 
quet Hamiltonian [33j . The explicit form of this operator 
is given by 



-v h in O 

H m (t) = --[(J z + iA?(t))OT(t) + h.c] — 
+2^[(J z+l A T (t)fd" 2 l (t) + h.c} 

-^[^ + (Ar(*)) 2 ]-f (K(m 2 - 



We consider here the notation 



6?(t) = E J i 

l — ~oo 

oc 

62 (t) = E & 



^-(2A™(<) + 1) 



^(AS*(*) + 1) 



tint 



iint 



(12) 

(13) 
(14) 



where Ji(z) is the /th-order Bessel function [37|, and 

A\ n (t) = -J y cosQ m (t) - J x sin0 m (i), (15) 
A5*(t) = J x cos6» m (t) - J y swO m (t). (16) 

The Hamiltonian Eq. ([T2|) can be written in the form 

00 

H m (t)= E exp (inOi)- (17) 

n— — 00 

In analogy with the standard RWA of quantum optics, we 
obtain an approximate Hamiltonian to describe the mth 
resonance by neglecting all the terms in H m (t) with os- 



cillatory time-dependence: H rn (t) 



The effective 



Hamiltonian governs the dynamics in the rotating 
frame. 

We next perform a bosonization procedure of via 
the Holstein-Primakoff representation Eq. ([3]). In the 
bosonized version, the effective Hamiltonian is written in 
terms of the bosonic mode b. To investigate the criti- 
cality in the system we introduce a complex macroscopic 
displacement of order for the bosonic operator as 
follows 



(18) 



c + aVN, 



where a = (Q + iP) (Q and P arc dimcnsionlcss param- 
eters) and c is a bosonic operator describing quantum 
fluctuations in the system. In the thermodynamic limit 
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N — > oo, wc perform a series expansion of the effective 
Hamiltonian in powers of yiV 

= ^ m) (c,c+) + VF ^ m) (c,ct) + NE^\Q,P), 

(19) 

where hp^^c, c 1 ) is a quadratic bosonic Hamiltonian (dc- 



Q 

pending on the choice of the macroscopic displacements 
Eq. (HI), h^< 

E { g\Q,P) is the lowest quasienergy (LQE). 



h\" (c, c^) contains linear bosonic terms, and 



The effective Hamiltonian for the m — case 



Instead of performing an abstract general theory for 
the effective Hamiltonian related to general m-photon 
resonances, we focus here on an illustrative particular 
case that contains the more relevant information, i.e., 
we consider the case m = 0. In this case, the effective 
Hamiltonian reads 



/i (0) - 
ri — 



iJy)J0 



^(2J 

Nn [ Jx 



Nn [Jx 



+ h.c 



-h^l + JD-^Jl- (20) 



2N 



N 



In the thermodynamic limit N — > oo, we expand the 
Holstein-Primakoff representation Eq.Q with respect to 
the complex macroscopic displacement Eq. (|18j) up to 
leading order N 



J x = NQ^l~\a\\ 



J V =NP\ 1 



J Z =N\ \a\ 



(21) 
(22) 

(23) 



This transformation describes a mapping from the co- 
ordinates (Q, P) onto the Bloch sphere, because the 



norm of the spin vector (J x , J y , J z ) has constant length 
J = N/2. From Eqs. and IJ22J it follows that 

\a\ 2 < 1. Furthermore, one can see from these rela- 
tions that all points (Q,P) with \a\ =1 correspond to 
J z = N/2. Therefore, all points of the boundary \a\ =1 
are mapped into the north pole of the Bloch sphere. For 
the interior points |a| 2 < 1, the transformation is bijec- 
tive, e.g., the origin (Q,P) = (0,0) is mapped into the 
south pole. 

By replacing Eqs. Ij2l ]) .(p2 ]) . and (J23J) into the effective 
Hamiltonian Eq. (|20|) we obtain the LQE for the m = 

case 



E ( g\Q,P) 



M 2 - \ ) 



2t| 



i-I«I 2 )p 2 



Jo 



a 



i-I«I 2 )p 2 



7 j(l-|a| 2 )g 2 . (24) 



r 



In the undriven case, a classification of the quantum 
phases of the LMG Hamiltonian is usually performed 
by studying the global minima of the energy landscape 
[15| . However, in this work we interpret the local minima 
of the LQE as metastable phases of the driven system. 
These metastable states are related to the phenomenon of 
parametric stabilization (32[. Since these are separated 
from the global minima by macroscopic displacements, 
we expect transitions to be suppressed, such that the 
corresponding values of the order parameters are observ- 
able. This possibility is reinforced when one recalls that 
Eq(Q,P) docs not have the same thermodynamic sig- 
nificance as the lowest actual energy. Therefore, in this 



work we consider the number of minima of the LQE land- 
scape as a criterion to establish a new phase diagram. 

FigH] (a) depicts the phase diagram as a function of 
7 y and jf , for fixed 7q/|/i| = 0.5, and FigfT] (b) as a 
function of 7q and 7 f , for fixed r y v /\h\ = 2. In these 
phase diagrams, we see many regions corresponding to 
different number of minima of the LQE landscape. Ad- 
ditionally, the phase diagrams show the appearance of 
many novel metastable phases, which are separated from 
each other by boundary lines, whose crossings correspond 
to noncquilibrium multicritical points. By crossing the 
line 7 y = —h in Figfl] the single global minimum at 
(Q, P) = (0, 0) splits into two macroscopically separated 
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FIG. 2: Lowest quasienergy landscape for the parameters 
]77| {h, 7o , 7i ,7 y ) = (-1, -1, 210, 2). In the undriven system, 
these parameters correspond to the symmetry-broken phase. 
The quantum evolution for N = 100 particles within one pe- 
riod is calculated when the system is initialized in a spin co- 
herent state. The green (light gray) line on the surface depicts 
the evolution of an initial wave packet centered at minima B 
and the red line shows the corresponding evolution for a wave 
packet initially centered at A. We consider the parameters 
Q/\h\ =40 and h/\h\ = -1 



global minima thus resembling the second order QPT 
known from the time-independent model. Interestingly, 
regions with even and odd number of minima are char- 
acterized by the existence of two and one global minima 
respectively. We can study the stability of the global 
minimum at (Q,P) — (0,0) analytically by computing 
the Jacobian-matrix at the origin of the LQE landscape 
and its eigenvalues 



Ax = -2(h + i»), 

A 2 = -2h - 2 7o * -(h + f) 



(25) 
(26) 



Both eigenvalues are positive in the region 7y < —h (the 
contours Ai = and A2 = are depicted as red (light 
gray) and black dashed lines respectively in FigO}. In 

the region 7 y > — h and jf < flJ (region be- 
tween the dashed curves in Fig[T](a)), Ai is negative and 
A2 is positive. Furthermore Ai and A2 are negative for 

7f > fl^J [jz^zT^j 7 and therefore, by crossing the curve 

A 2 = 0, the saddle point at the origin becomes a local 
maximum. Consequently, in the region 7 y > — h the sin- 
gle central minimum splits up in two global minima. The 
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FIG. 3: Quantum evolution of the observables for finite size 
N = 100 and parameters 7o , 7?, 7 H ) = (-1, -1,210,2). 

The expectation values {J x )/j, (J y ) / j and {Jz)/j are depicted 
by the orange, cyan and magenta curves respectively, (a) De- 
picts the quantum evolution within one period and (b) depicts 
the stroboscopic dynamics when the system is initialized in 
a wave packet centered at the minimum A depicted in fig- 
ure [2] Correspondingly, (c) and (d) depict the inter-period 
and stroboscopic evolution, respectively, when the system is 
initialized in a wave packet centered at the minimum B de- 
picted in figure [2] We consider the parameters Q/\h\ = 40, 
h/\h\ = -1. 



phase diagram depicted in FigfT] (b) is characterized by 
Ai < 0, and therefore, the dashed line corresponds to the 
boundary between the regions below and above the level 
curve A 2 = 0, where the origin is a saddle point and a 
maximum respectively. 

An example for the LQE is given in Fig. [3J where 
the parameters are so chosen, that the undriven sys- 
tem (7* = 0) is in the symmetry-broken phase. As the 
LQE corresponds to the ground state energy in the un- 
driven case, it exhibits two global minima corresponding 
to macroscopically separated states degenerate in energy. 
In the driven system, apart from the two global minima 
characteristic of the undriven system (in Fig. [5] denoted 
with A), new characteristics of the LQE landscape ap- 
pear, e.g. two local minima (denoted with B). These new 
local minima can be interpreted to be novel metastablc 
states, which are strongly related to the quantum evolu- 
tion of the system as we describe in the next section. 



B. Quantum evolution 

In this section we investigate the quantum evolution 
when the system is initially prepared in a spin coher- 
ent state |38l | centered at a local minimum of the LQE 
surface. Spin coherent states have minimum uncertainty 
and are the closest ones to a classical angular momentum 
state. To describe geometrically the quantum evolution, 
we parametrize the Bloch sphere using spherical coordi- 
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nates and express the angular momentum components in 
terms of the azimuthal (<fi) and polar (0) angles as follows: 



N 

J x = — sin 9 cos <j), 
N 

J y = — sin 9 sin <p, 



Jz 



N 



cos 9. 



(27) 
(28) 
(29) 



By replacing this set of equations into Eqs. ([2T|) -([23 | we 
find a relation between the (9, <j>) and (Q, P) coordinate 
systems 



arccos 



arccos [2 (Q 
2Q 



P 2 ) 



1] 



sin (9) 



(30) 
(31) 



With the angular coordinates <f> and 9, we can repre- 
sent the spin coherent state \</>, 6} by using \<f>,8) = 
(1 + |r| 2 )-^ exp[r J+] \j,-j) with r = e^tanf . This 
procedure describes a mapping from a point of the LQE 
surface onto the set of spin coherent states. 

To study the quantum evolution we consider a sys- 
tem consisting of N = 100 particles. For a finite size 
N, the numerical problem consists in the solution of 
2N + 1 coupled ordinary differential equations. After 
the numerical integration of the Schrodingcr equation, 
we construct the evolution operator U(t, 0), which allows 
us to calculate the state of the system at any time t > 0: 
\^,t) = U(t,0) \<t>,9), providing that we prepare initially 
the system in a spin coherent state |^,0) = \(j>,&). Wc 
now proceed to calculate the normalized expectation val- 
ues of the spin components in the state Y&,t). Fig. [3] 
(a) depicts the continuous time evolution within one pe- 
riod of the driving T = 2tt/Q of an initial wave packet 
centered at the minimum A in Fig. [2] Similarly, Fig. 
[3] (c) shows the inter-period dynamics when the initial 
wave packet is centered at the minimum B in Fig. [5] 
To obtain a better geometrical picture of the quantum 
evolution we represent the mean values of the angular 
momentum components by means of the coordinates Q 
and P. Therefore, we project the expectation values onto 
the Bloch sphere by calculating the angles <j> and 9 and 
then solving Eqs. $U$) and for Q and P. The re- 
sult is shown in Fig. [5] for N = 100 particles, where 
the time evolution for initial wave packets centered in 
minima A and B is depicted by the red and green (light 
gray) curves, respectively. The trajectory initialized in 
A is strongly trapped within the minimum, whereas the 
other trajectory exhibits higher oscillations around the 
initial state. For finite size A^ ^> 1 the trajectories take 
place approximately over the surface of the LQE land- 
scape, i.e., the mean value of the spin evolves in average 
along points with equal values of E G . This behavior is 
connected to the fact that the eigenvalues of the effec- 
tive Hamiltonian — from which the LQE is derived, corre- 
spond to the quasiencrgics of the system and the average 
value of the quasienergy is conserved in a time-periodic 



system. The state in the laboratory frame \^f,t) and the 
state in the rotating frame \^ m ,t) — in which the LQE 
is derived, are connected via 



|*,*) =U m {t)\* m ,t) 



U m (t)e 



-ill 



£(m) f 



(32) 
(33) 



where U rot (t, 0) 



-ih). 



denote the propagator in 



the rotating frame and |$ m ,0) = 1^,0) as a conse- 
quence of Eq. (fTTj) . Furthermore, the propagators in the 
laboratory and rotating frame are related as U (t, 0) w 



U m (t)t 



Additionally, there is a very interesting 



relation between the stroboscopic quantum evolution and 
the parity operator Eq. ^ that can be established using 
the relation 



Um (tr ) 



fl exp 



ttN 



(34) 



for t r — 2 j^ L with integer r. As a consequence of this, 
for m = we have U m (t r ) = /, which implies that the 
states in the laboratory frame and in the rotating frame 
are identical for these times. Thus, the stroboscopic time 
evolution is governed entirely by . The stroboscopic 
long-time evolution of the observables is displayed in Fig. 
[3fb). Here, the quantum evolution is strongly trapped in 
the neighborhood of the minimum A in Fig. [5] A similar 
situation occurs in Fig. Eljd) for a wave packet initially 
centered at the minimum B in Fig. [2j By calculating 
the quantum evolution with different system sizes, one 
finds that the quantum fluctuations of the trajectories 
decrease with number of particles, and therefore, in the 
thermodynamic limit, the time evolution has to be con- 
strained to the LQE landscape. Because the time evo- 
lution is connected with the geometrical features of the 
LQE landscape, it is justified to use it as a background 
to define the existence of new metastable phases. 



III. CONCLUSIONS 

We have investigated the driving-induced QPT in the 
LMG model. We show that under the effect of an exter- 
nal driving the system exhibits a multistable character 
with no analogue in equilibrium systems. In particular, 
the novel quantum phases correspond to local minima 
of the LQE landscape. To understand the nature of the 
noncquilibrium metastable states, we study the quantum 
evolution for finite size N when the system is initially pre- 
pared in a coherent state centered at one particular local 
minimum of the LQE surface. We find that the system 
is dynamically trapped in the neighborhood of the initial 
state, and the quantum evolution is related to the local 
curvature of the chosen minima, i.e., for higher curvature 
the trapping effect is stronger. Our approach opens a new 
window in the understanding of driving-induced critical- 
ity, in particular, the effective Hamiltonian Eq. (|20[) con- 
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tains effective interactions which are absent in equilib- 
rium, in this sense, the effective Hamiltonian can be con- 
sidered as a quantum simulator. Surprisingly, despite of 
the close relation between the LMG model and the Dickc 
model, in contrast to our description of the nonequilib- 
rium QPT in the Dicke model [13j , the order of the phase 
transitions in the LMG model does not change under the 
effect of driving. In the case of undriven QPTs, the crit- 
ical behavior is related to nonanaliticities of the ground 
state. Under the effect of an external driving, however, 
the system will experience transitions to excited states, 
even if it is initially prepared in the ground state. In 



this sense, the global minima as well as the local minima 
of the LQE landscape play the role of driving-induced 
states of matter. 
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